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Abstract: A pattern matching based tracking algorithm, named MdcPatRec, is used for the reconstruction of charged tracks 
in the drift chamber of the BESIII detector. This paper addresses the shortage of segment finding in MdcPatRec algorithm. 
An extended segment construction scheme and the corresponding pattern dictionary are presented. Evaluation with Monte- 
Carlo and experimental data show that the new method can achieve higher efficiency for low transverse momentum tracks. 
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1 Introduction 

The Beijing Electron Positron Collider II (BEPCII) H 
is a double-ring multi-bunch e + e~ collider which is oper- 
ated in the collision energy range from 2 GeV to 4.6 GeV. 
The Beijing Spectrometer III (BESIII) [1] is a general- 
purpose detector with 93% coverage of full solid angle. 
BEPCII and BESIII are powerful facilities for the study 
of charmonium physics, D-physics, hadron spectroscopy 
and r physics Q2|] . From the interaction area of the e + e~ 
beams to the outside, there are 5 apparatuses: a drift 
chamber (MDC) y|], which includes 43 layers of drift 
cells and is used for charged tracking; a time-of-flight 
counter, which is composed of 2 layers of scintillator bars 
for measuring the flight time of penetrating charged par- 
ticles; an electromagnetic calorimeter, which comprises 
6240 CsI(Tl) crystals for the measurement of photon en- 
ergy; a super-conducting solenoid magnet, which provides 
a 1 Tesla magnetic field and a muon counter, which in- 
cludes 1000 m 2 resistive plate chambers, used for identi- 
fying muon. 

MDC is a carbon fiber cylinder flowed with a mixture 
of He and C 3 H 8 (He/C 3 H g = 60/40) as working gas, and 
with signal and field wires strung between the two end- 



caps. The signal wires are operated with a high positive 
voltage, while the field wires are at ground voltage. There 
are 6796 signal wires in total; each of them is surrounded 
by 8 or 9 field wires. The signal wires are grouped into 43 
layers, and these layers are divided into 1 1 super-layers, 
with the labels SL- J (from the inside to the outside, J = 
1,2, • • • , 11.). The SL-1 1 consists of only 3 signal layers, 
while each of the other 10 super-layers includes 4 signal 
layers. In some super-layers, the signal wires are parallel 
to the z axis of the MDC, so these super-layers are named 
axial super-layers. In the other super-layers, the signal 
wires have a stereo angle with respect to the z axis of 
the MDC, and these super-layers are named stereo super- 
layers. Axial (A) and stereo (U V, U means the stereo 
angle is negative, V means the stereo angle is positive) 
super-layers follow the order UVAAAUVUVAA fl. The 
small space around a signal wire and surrounded by the 
neighbouring field wires is called a drift cell. The sec- 
tions of these drift cells are nearly square in the plane 
perpendicular to the beam axis, with the size of approx- 
imately 12 x 12 mm 2 in the innermost two super-layers 
and 16 x 16 mm 2 in the other super-layers. The signal 
wires collect the avalanche electrons produced by ioniza- 
tion along the particle trajectories. Electric signals on fired 
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wires are called hits and are used to reconstruct trajecto- 
ries of charged particles. 

One basic offline C++ software algorithm for MDC 
track reconstruction at BESIII is a pattern matching based 
tracking package (MdcPatRec) |4], which extracts track 
segments from a set of hits within a group of neighbor- 
ing cells and links the segments into tracks. For charged 
particles with high transverse momentum (p t ), the track- 
ing efficiency of MdcPatRec is satisfied for physics anal- 
ysis. However, in low p t range, the efficiency still has a 
large space to improve. In this paper we propose an ex- 
tended segment construction scheme and its pattern dic- 
tionary which is justified to be helpful in the improvement 
of low p t tracking efficiency. 

2 MdcPatRec algorithm 

To introduce the procedure of MdcPatRec, we list the 
important glossaries first: 

(1) Segment. In a super-layer, a series of hits of a 
charged particle are defined as a segment. The arc in Fig. 
Q] shows a track. According to the segment definition, the 
4 hits at the 4 fired wires compose a segment. 
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a segment and an 8-wire group 
Fig. 1. A schematic plot for illustrating a segment 
and an 8-wire group. The 4 layers of circles denote 
the signal wires in one super-layer. The arc shows a 
track, and the black solid circles are the fired wires 
of the particle. The 8 adjacent signal wires marked 
with a cross compose an 8-wire group used to match 
the segment pattern in segment finding. 

(2) Segment pattern and segment pattern dictionary. 
The fired wire distribution of a segment in a super-layer 
is defined as a segment pattern. In present MdcPatRec, 
segment finder uses the hits in the 2nd layer of a super- 
layer as seeds to search for segments in every group of 8 
signal wires, the relative position of these signal wires is 
graphically represented by the 8 wires marked by a cross 
in Fig. [T] A segment pattern is expressed by 8 bits of an 
integer in C++ program: the fired wires are indicated by 
"1", while the non-fired wires are indicated by "0". For 
instance, the segment pattern illustrated in Fig. [TJis ex- 
pressed with 01001101 in the binary number (77 in the 
decimal number). Segment patterns are collected to build 
a segment pattern dictionary. 



The procedure of MdcPatRec includes the following 
steps: first to find segments in all super-layers by pat- 
tern matching; then to assemble segments in axial super- 
layers into circular tracks and apply a circle fit. Second, 
to add segments in stereo super-layers to circle tracks to 
constitute helix tracks and apply helix fit. After helix fit, 
tracks are reconstructed, and they are stored for Kalman 
fit, tracks extrapolation, physics analysis and so on. 

3 The shortage of 8-wire segment pattern in 
MdcPatRec 

The tracking efficiency of MdcPatRec is satisfied for 
non-curled tracks in the drift chamber. It loses some effi- 
ciency for low p t tracks due to the reduced efficiency for 
segment finding within a super-layer. Azimuth coverage 
angle of segment groups does not meet the requirements 
for low p t track segment finding. 

Drift cell size limits the azimuth coverage of segment 
pattern dictionary. We compare the azimuth coverage of 
an 8-wire segment group with the azimuth coverage of 
tracks in a super-layer to test whether the segment group 
is large enough to reconstruct all possible segments. We 
define a variable Ar/> to describe the coverage of a segment 
group along cf> direction, the definition of A</> is described 
by Eq. © 



->L~L(L-R) ' 



>J U-R(JJ-L) |) 



(1) 



where <fi L - L is the 4> angle of the drift cell at the low left 
corner of an 8-wire group, as displayed in Fig. [T] simi- 
larly, L— R means the low right corner, U — L means the 
upper left corner, and U — R means the upper right cor- 
ner. The A(f> values of 8-wire segment pattern dictionary 
are listed in Table [TJ We also define a variable, 5(f>, to 
describe the azimuth coverage of a track in every super- 
layers along 4> direction. In the magnetic field, 5<fi of a 
track with certain p t within a super-layer can be calcu- 
lated by Eq. ©. 
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In Eq. (0, r in and r out are the inner radius and outer 
radius of a super-layer respectively; r v is the radius of a 
track with certain p u which can be calculated by Eq. (f3]) 

Pt 



B z -q 



(3) 



where B z is the magnetic field intensity along the beam 
axis and q is the electric charge of the particle; all vari- 
ables are expressed in the international system of units. 
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The calculated 5(f) are listed in Table [2] Comparing these 
5(f) with A(f> listed in TableQ] we find that the 8-wire seg- 
ment pattern dictionary is not large enough to cover all 
segments of some low p t charged particles. 

The wire groups used to search for segments in the 
MDC are not all symmetric along <f> direction. If the 
amounts of wires of two adjacent layers are the same, the 
two layers are staggered exactly by a half drift cell along 
<t> direction, as shown in Fig. Q] Otherwise, the two layers 
cannot be placed in the exact staggered form. As a result, 
the symmetric arrangement of the wire groups in a super- 
layer is determined by the amounts of wires in the four 
layers. From SL-1 to SL-5, the amounts of wires in each 
layer are not the same, and some 8-wire wire groups in 
these super-layers are asymmetric. Figure shows two 
asymmetric 8-wire wire groups in SL-1. As shown in 
these plots, the azimuth coverage of an asymmetric 8-wire 
wire group is smaller than that of a symmetric one. Table 
Q]lists the minimum wire group coverages of 8-wire seg- 

Table 1. The azimuth coverages, i.e. A(f> (degree), of 8-wire and 14- wire segment pattern dictionary in each super-layer. 
The 'min' value is the minimum coverage in a super-layer. The wire groups are not all symmetric in SL-1 to SL-5, so 
the minimum A(f> is smaller than the symmetric wire groups. In SL-6 to SL-11, all wire groups are symmetric, so the 
minimum A<j) is equal to that of the symmetric wire groups. 
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Table 2. The track azimuth coverage, i.e. S(f> (degree), of charged particles with different pt in every super-layer. If the 
pt of a particle is smaller than 120 MeV/c, the particle curls in the MDC. If the 8(f> is larger than the minimum coverage of 
8-wire segment pattern dictionary in this super-layer, it is marked with an underline. If the 5(f) is larger than the minimum 
coverage of 14-wire segment pattern dictionary in this super-layer, it is marked with an overline. 
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ment pattern dictionary in every super-layer; as expected, 
the asymmetry of wire groups in SL-1 to SL-5 reduces the 
coverage of segment pattern. In addition, the 8-wire seg- 
ment pattern dictionary is designed based on the hypoth- 
esis of symmetric wire groups, and the dictionary doesn't 
suit the asymmetric wire groups. 




-3-2-10123 -3-2-10123 



x (cm) x (cm) 

Fig. 2. The relative position of signal wires of 2 
asymmetric 8-wire groups in SL-1. The origin of 
coordinate in each plot is set at the position of the 
signal wire in the second layer. 



4 The extended segment pattern dictionary 

In Sec. 3, we explain 8-wire segment pattern dictio- 
nary which is not optimal in the azimuth coverage for 
tracks with low p t . In addition, the asymmetry of the wire 
groups aggravates this problem. Hence, it is meaningful 
to broaden the 8-wire groups to larger wire groups and de- 
sign an extended segment pattern dictionary for segment 



finding. 

An easy way to broaden the coverage of the 8-wire 
segment pattern dictionary is to add a signal wire at each 
side of an 8-wire group on the 1st, the 3rd and the 4th 
layer, and an 8-wire group is extended to a 14-wire group 
as illustrated in Fig. [3] Comparing the azimuth coverage 
of particle tracks in Table |2] with the coverage of 14-wire 
segment pattern dictionary in Table [T] we can find that the 
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14- wire segment pattern dictionary can cover all segments 
of different tracks except the segments in the outermost 
super-layer, which a track can reach. This type of seg- 
ments has very large azimuth coverage angle and cannot 
be reconstructed in an easy way. 

A new segment pattern dictionary is built based on the 
14-wire groups. Because the least squares fitting is ap- 
plied to the hits of a segment pattern in segment finding, 
a segment pattern requires 3 or more hits. Therefore we 
build the dictionary with 3-hit and 4-hit segment patterns. 
A 4-hit segment pattern has a hit in every layer, and a 3- 
hit pattern has 3 hits in three different layers. For the 4-hit 
segment patterns, there are 80 (4 x 1 x 4 x 5) combinations, 
but some combinations are obviously impossible and we 
remove them, which are determined by the relative posi- 
tion of the signal wires in a 14-wire group. From SL-6 
to SL- 11, all wire groups are symmetric, and we build a 
unique segment pattern dictionary for them. There is no 
significant asymmetry from SL-2 to SL-5 and we find that 
the dictionary based on symmetric wire groups is appli- 
cable. In SL-1, the asymmetry of wire groups is more 
serious than any other super-layers, and we build a special 
segment pattern dictionary. Finally the segment pattern 
dictionary for super-layers from SL-2 to SL-11 includes 
30 4-hit segment patterns and 56 3-hit segment patterns, 
and the segment pattern dictionary for SL-1 includes 44 
4-hit segment patterns and 72 3-hit segment patterns. Fig. 
|4]shows 6 4-hit segment patterns. In C++ program, a 14- 
wire segment pattern is expressed by 14 bits in an integer 
variable: the fired wires are marked by "1", while the other 



signal wires are marked by "0". Table |3]lists the values of 
decimal numbers corresponding to the segment patterns 
drawn in Fig. |4] 

O O ^ 3 ^ 2 ^ 11 ^ ^ O O 4th layer 
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OOOO^OOOO 2nd layer 
O O^^^^O O 1st layer 

a 14-wire group 

Fig. 3. A 14-wire group schematic plot. The 4 
layers of circles denote the signal wires in a super- 
layer. The 14 adjacent signal wires marked with a 
cross make up the 14-wire group. 
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Fig. 4. Some 4-hit 14-wire segment patterns. 



Table 3. The expression of the 4-hit segment patterns listed in Fig. [4] 
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5 The 14-wire segment pattern dictionary 
performance 

With single charged particle Monte-Carlo (MC) data 
sample, a multi-prong physics MC data sample and exper- 
imental data sample, the performance of 14-wire segment 
pattern dictionary is tested. 

5.1 The segment efficiency 

The segment efficiency e seg in a super-layer is defined 
as Eq. © 

jyrcc 

= scg (4) 



where iV s r c e g is the number of segments that have been 
found in a super-layer, and is the number of tracks 
passing through the super-layer in MC data sample. Fig- 
ure[5]shows the segment efficiency of e~ with different p t , 
as shown in these plots, the 14-wire segment pattern dic- 
tionary can achieve higher segment efficiency at the super- 
layer which track curled back. 

5.2 The tracking efficiency 

We generate 400,000 single e~ MC data sample to test 
the improvement of tracking efficiency. In the MC data 
sample, the momentum of e~ has a uniform distribution 
from 0.05 to 1 GeV/c, and the angular distribution of the 
er is uniform in full solid angle. Tracking efficiency e tr k 
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is defined as Eq. (O 



£trk 



(5) 



where N™£ is the number of events that have recon- 
structed a charged track, and A^ 4 k c is the number of events 
that have a charged track in MC data sample. Figure [6] 
shows the tracking efficiency for e~. As shown in these 
plots, the tracking efficiency of 14- wire segment pattern 
dictionary is higher than that of the 8-wire segment pat- 
tern dictionary. 
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Fig. 5. The segment efficiency of e~ with different 
pt in each super-layer. Symbol x and + are used to 
mark the segment efficiencies of 14- wire and 8-wire 
segment pattern dictionaries respectively. 
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Fig. 6. The tracking efficiency of e~. 
5.3 The transverse momentum resolution 

Using the same MC data sample we generate for test- 
ing the tracking efficiency, we compare the transverse mo- 
mentum resolution between the extended pattern dictio- 



nary and the old pattern dictionary. The transverse mo- 
mentum resolution 5p t is an index of track reconstruction 
quality. S pt is expressed by Eq. (O 



5p t = PT 



-Pt 



(6) 



where p\ cc is the p t of a reconstructed track, while pf c is 
the p t of this charged particle in MC simulation. Fig. [7] 
shows the 5p t distribution of e~ (p t less than 130MeV/c) 
reconstructed by the 14-wire and 8-wire segment pattern 
dictionaries respectively. By fit, the 5p t resolution is im- 
proved from about 0.94 MeV/c to about 0.89 MeV/c. 
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Fig. 7. The Spt distribution of e~. The two his- 
tograms are normalized. 
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Fig. 8. The hits number distribution of e . 

5.4 The noise resistant ability 

With the same MC events that we use to test the track- 
ing efficiency, we test the noise resistant ability of the ex- 
tended pattern dictionary. The track hits number is the 
amount of hits used in the reconstructed tracks, and it de- 
termines the quality of track reconstruction. The track 
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ability of tracking algorithm. The definition of R n 
expressed by Eq.(|7]i 



R„ 



iV hit 
iv hit 



(7) 



where N£° t mc is the number of fired wires used in track 
reconstruction but they are not fired by the particle; jVJ°* 
is the number of signal wires used in track reconstruction. 
Figure [8] and [9] display the number of hits on track and 
^noisc of e~ with p t less than 130MeV/c. In these plots, 
the average N^ft is increased from 33.1 to 34.9, while the 
average i? no i se is decreased from 0.43% to 0.42%. The 
change of track noise level is very little. 
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Fig. 9. The R no i BC distribution of e 
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Fig. 10. The reconstructed signals of r\ from MC 
decay sample ip(2S) —> "iXC.J — >• r )riK + K~ with 
the requirement p t of e~ less than 130MeV/c. 

5.5 Application in physics analysis 

We generate a physics process, i/j(2S) —> jXcj 
^fK + K~rj — > jK + K~e + e~, to test the performance of 
the new segment pattern dictionary. Figure [10] is the in- 
variant mass distribution of rj. As shown in the plot, the 



new segment pattern dictionary can keep more physical 
events than the old segment pattern dictionary. The total 
reconstruction efficiency is increased by about 4%; if we 
only count the events with p t of e~ less than 130 MeV/c, 
the efficiency is increased by about 20%. 

5.6 Application in experimental data 

We test the new segment pattern dictionary with ex- 
perimental data. Figure QT| shows the number of recon- 
structed tracks with respect to p t for a sample of ^(3770) 
experimental data. The amount of reconstructed tracks 
with p t less than 130MeV/c is increased by about 15%. 
The number of charged tracks of a ^(3770) event is usu- 
ally not less than 4. The differences in the number of 
charged tracks and in the track angular distribution be- 
tween ^(3770) events and single e~ MC events make the 
difference between the curves of reconstructed tracks in 
Fig. QT|and efficiency curves in Fig. [6] 
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Fig. 11. The reconstructed tracks numbers from 
t/>(3770) experimental data. 

5.7 CPU-time consumption 

Within BESIII offline software system Jl], We test the 
CPU-time increase of MdcPatRec using the 14-wire seg- 
ment pattern dictionary. Compared with using the 8-wire 
segment pattern dictionary, the CPU-time of MdcPatRec 
in segment finder is increased by about 2 ms for an event 
on the average; the CPU-time in total MdcPatRec algo- 
rithm is increased by about 30%; the CPU-time of the total 
event reconstruction is increased by about 10%. 

6 Summary and discussion 

In this paper, we study the issue of low tracking ef- 
ficiency of the MdcPatRec algorithm in the low p t range 
and understand one source of the problem as the short- 
age of 8-wire segment group: its limit of azimuth angle 
coverage and wire group asymmetry. We propose a 14- 
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wire group and its segment pattern dictionary as a rem- 
edy to the defect. The performance shows that the new 
wire group construction scheme and its segment pattern 
dictionary can improve the tracking efficiency for the low 



p t tracks. However, the CPU-time consumption is in- 
creased, so we will do more work to optimize the algo- 
rithm. 
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